#ifndef LINEAREQNSOLVER_H_
#define LINEAREQNSOLVER_H_

#include<iostream>
#include<vector>
#include<Eigen/Dense>

// 三对角矩阵系数(l下次对角线, u上次对角线, d主对角线)
std::vector<double> solu(const std::vector<double> &_l,
                        const std::vector<double> &_u,
                        const std::vector<double> &_d,
                        const std::vector<double> &_b);

// 一般矩阵系数
std::vector<double> solu(const std::vector<std::vector<double>> &_A, 
                        const std::vector<double> &_b);

void print_vector(const std::vector<double> &x);

void print_matrix(const std::vector<std::vector<double>> &A);

#endif